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Resumen 


El estudio fue realizado en el periodo de lluvias de 2006 en tres sitios del semiárido 
centro de México: “El Carmen”, en el estado de Guanajuato; “Amazcala” y 
“Cadereyta”, en el estado de Querétaro; ello, con el objetivo de estudiar el 
comportamiento de las pérdidas por intercepción, E, de dos especies dominantes, 
mezquite (Prosopis laevigata) y huizache (Acacia farnesiana). El experimento se basó 
en la metodología de Guevara-Escobar et al. (2007). La modelación de las E, por 
tormenta fue construida con dos modelos de base física, Rutter, Gash, y un tercero 
del tipo lineal múltiple (RLM). En cada sitio fueron seleccionados dos árboles de 
ambas especies. La precipitación en 2006 fue de 770 mm en El Carmen y de 732 
mm en Cadereyta. En Amazcala, en el periodo agosto-noviembre, de 451 mm. 
En proporción a la precipitación, las E, por evento resultaron de 15, 20 y 27% en 
huizache, y de 27, 22 y 17% en mezquite, para Cadereyta, El Carmen y Amazcala, 
respectivamente. Estadísticamente, las eficiencias de modelos de Gash, Rutter y 
RLM correspondientes a los tres sitios resultaron de -0.1<E<0.34, -4.41<E<0.8 y 
-1.19<E<0.94. Las diferencias en porcentaje entre las observaciones y los modelos 
para el periodo de estudio fueron de <120%, siendo el modelo RLM el más eficiente 
(0%). La prueba ANOVA entre árboles y sitios no mostró diferencias estadística- 
mente significativas, por lo cual fue generada una correlación para estimar las E, en 
el contexto regional, de donde se obtuvo un promedio del 20% de intercepción. 


Palabras clave: pérdidas por intercepción, evaporación, percolación, flujo 
cortical, balance hídrico, ecosistemas semiáridos y modelación. 


Introducción 


La vegetación es el primer obstáculo en 
interceptar la precipitación y modificar su 
trayectoria antes de llegar al suelo. Desde 
el punto de vista hidrológico, es un estrato 
filtrante de la precipitación que afecta la 
repartición de la evaporación, la infiltración 
y los escurrimientos superficiales; por tanto, 
es de gran importancia en la distribución 
espacial del balance hídrico (Rutter et al., 
1971; Schellekens et al., 1999; Zeng et al., 2000; 


Van-Dijik y Bruijnzeel, 2001). El proceso de 
intercepción o pérdidas por intercepción es un 
concepto donde se involucra cualquier tipo de 
cubierta sobre el suelo, sea natural o artificial. 
Integra la precipitación retenida y reintegrada 
a la atmósfera en forma de vapor en todos los 
estratos vegetales de un ecosistema (arbóreo, 
arbustivo, herbáceo), incluso el agua retenida 
en la capa superficial del suelo formada con 
residuos vegetales, producto de la senescencia 
de la vegetación, ampliamente conocido como 
mulch (Groen y Savenije, 2006). 
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En general, las pérdidas son controladas por 
la arquitectura de la vegetación y las variables 
(Crockford y Richardson, 
investigaciones 


climatológicas 
2000). Las 
bosques templados y tropicales se han hecho 
con diversas técnicas. 


efectuadas en 


Habitualmente, los 
estudios reportan las pérdidas en proporción 
de la precipitación anual. Las pérdidas han 
sido ampliamente estudiadas en bosques 
fragmentados y en bosques de reserva ecológica 
de climas templados y tropicales (Rodrigo y 
Ávila, 2001; Loescher et al., 2002; Raat et al., 
2002). Sin embargo, aún falta conocer más sobre 
el proceso de las pérdidas por intercepción 
en ecosistemas de ambientes semiáridos. De 
la misma forma, ampliar el conocimiento en 
árboles aislados y zonas de transición entre 
dos ecosistemas distintos o ecotonos, donde la 
interacción entre las fronteras o bordes influye 
en la distribución espacial de la intercepción. 

Los trabajos experimentales y de 
laboratorio llevados a cabo para estudiar el 
proceso de intercepción en árboles aislados, 
bosques naturales y urbanos, han permitido 
desarrollar diversos modelos para simular 
el proceso de intercepción (Horton, 1919; 
Rutter et al., 1971; Gash, 1979; Liu, 1997; Xiao 
et al., 2000). No obstante, no hay modelos 
específicos para ambientes semiáridos, por 
lo que comúnmente suelen adaptarse los 
modelos existentes a este tipo de ambientes 
(Navar y Bryan 1994; Navar et al., 1999), 
principalmente porque la vegetación es 
poco valorizada y porque su composición 
estructural difiere de la de los bosques 
tropicales, templados o fríos. Es decir, la 
vegetación es dispersa, crece aislada o 
agregada, con especies dominantes en forma 
de islas, separadas entre sí por áreas de suelo 
desnudo o vegetación efímera; la interacción 
entre sus componentes es compleja y poco 
conocida su participación en la repartición de 
los componentes del balance hídrico. 

El objetivo del presente estudio fue evaluar 
el proceso de pérdidas por intercepción, 
mediante la cuantificación y modelación, 
basándose en un análisis por eventos en dos 


especies dominantes de la región semiárida del 
centro de México: huizache (Acacia farnesiana 
(L.) Willd.) y mezquite (Prosopis laevigata Humb. 
et Bonpl. ex Willd), con el supuesto de que las 
pérdidas a escala regional pueden expresarse 
por el comportamiento de árboles aislados de 
tres sitios ubicados dentro la región semiárida 
del centro de México. 


Teoría 
Intercepción 


La evaporación total en un ecosistema (ET) o 
cuenca la integran: 


ET =E¿+Eyy+ Ep +E¡+Ey, (1 


donde E, es la evaporación del suelo 
desnudo; E,,, la evaporación de los cuerpos 
de agua; En la evaporación por transpiración 
de la vegetación; E, la evaporación de la 
precipitación retenida por la vegetación, y E,, es 
la evaporación del agua retenida por el mulch; 
estas dos últimas ampliamente conocida como 
pérdidas por intercepción (Groen y Savenije, 
2006). Conceptualmente, la evaporación por 
las pérdidas por intercepción contiene dos 
componentes básicos: la evaporación a partir 
del inicio de la lluvia hasta la saturación de la 
vegetación y la evaporación desde la saturación 
hasta la terminación de la lluvia (Horton, 
1919): 


E t 
E¡=f Edt+ [Edt+S (2 
0 E 


donde tf es el tiempo de duración de la 
lluvia; f”, el tiempo para lograr la saturación 
de la vegetación, y S es la capacidad de 
almacenamiento del follaje o dosel. Las 
pérdidas por intercepción E, corresponden 
a la diferencia entre la precipitación total o 
incidente (P,), y la precipitación que llega al 
suelo o precipitación neta (P.). 


E¡=P5-P, (3) 
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La precipitación incidente, P,, sobre la 
cobertura vegetal produce dos tipos de flujo: 
flujo foliar o drenaje (hojas) y flujo cortical 
(tronco+ramas). Una parte de la lluvia es 
retenida y otra llega al suelo. Por lo tanto, 
la precipitación neta está constituida por la 
fracción de la precipitación que pasa libremente 
a través de la vegetación o follaje más la porción 
de la precipitación drenada después de ser 
retenida temporalmente. Así, las E, también 
pueden expresarse como: 


E¡=P¿-(T+ D+Sp) (4) 


siendo T, D y $, el flujo de percolación libre, 
drenaje foliar y flujo cortical, respectivamente. 
La integración de T y D suele denominarse flujo 
de percolación o traslocación, T,. 


Modelos 
Rutter et al. (1971 y 1975) 


En principio, el modelo de Rutter contiene cuatro 
parámetros asociados con la arquitectura de la 
vegetación, donde la cantidad de agua retenida 
temporalmente por la vegetación o el follaje es 
determinada por la percolación, la intensidad 
de la lluvia, la evaporación y el drenaje, de tal 
manera que: 


El p)R¡-E,-D para C>S (5) 


E U=p)R (0/5) E, =D paraC<S (6) 


donde C es la capacidad de almacenamiento 
de la vegetación o lámina de precipitación 
retenida por el follaje; R, la intensidad de la 
lluvia; S, el agua necesaria para humedecer 
toda la superficie de la copa o follaje; D, el 
drenaje del follaje; p, la percolación libre, y 
E, es la evaporación potencial estimada por 
la ecuación de Penman-Monteith (Monteith, 
1965). 


El drenaje del follaje es explicado mediante 
una función exponencial: 


D= D, e b(C=S) (7) 


siendo D, el drenaje cuando la copa está saturada 
C=S y b es un parámetro empírico. 


Gash (1979) 


Gash considera la lluvia representada por una 
serie de eventos discretos, suficientemente 
separados como para que la copa y el tronco de 
la vegetación puedan considerase secos entre 
cada evento. Con base en la expresión elemental 
de las pérdidas (ecuación (2)) por intercepción 
en proporción lineal con la precipitación tanto 
para eventos como para una serie de eventos se 
tiene: 


E, = ar, +b (8) 


Gash describe la ecuación (8) mediante las 
relaciones siguientes: 


b= 5+ [a Jo.) 


Al igual que Rutter et al. (1971 y 1975), Gash 
supone la capacidad de almacenamiento (S) 


Al] tm 


(9) 


en términos de la estructura de la vegetación 
y la define como lámina de lluvia retenida en 
la copa, en condición de cero evaporación, una 
vez concluida la lluvia y el drenaje. La lluvia 
drenada por el tronco y la drenada por las ramas 
la denomina p,. La percolación libre también la 
identifica con el término p. El agua almacenada 
en los troncos antes de iniciar el flujo cortical 
la llama capacidad de almacenamiento del 
tronco (S). Cuando considera la cobertura 
de la vegetación, a menudo asume p como el 
complemento de la fracción de la cobertura 
de la copa o la vegetación (c). Para una serie 
de eventos o un evento individual, el modelo 
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emplea la media de la evaporación (E) y la 
intensidad (R). 


Materiales y métodos 
Área de estudio 


La zona de estudio está situada en la región 
semiárida del centro del México; cubre una 
parte del corredor continental, constituido por 
la porción noreste del estado de Guanajuato y 
parte del centro-sur de Querétaro. La primera 
pertenece a la cuenca del Lerma-Santiago y 
la segunda, a la cuenca del Pánuco. El clima 
de la región es del tipo templado-semiárido, 
con un periodo de lluvias de verano con una 
lámina media anual entre 400 y 800 mm, y 
una temperatura media anual de 25 *C. La 
vegetación de la región la constituye el matorral 
xerófilo, representada en buena medida 
por pastos y una distribución dispersa de 
mezquite (Prosopis laevigata) y huizache (Acacia 
farnesiana). Los sitios de estudio se ubicaron en 
las microcuencas de El Carmen, Amazcala y 
Cadereyta, las dos primeras pertenecientes a la 


cuenca del río Laja, en los estados de Guanajuato 
y Querétaro, y la tercera situada en la parte alta 
de la cuenca del río Moctezuma, Querétaro 
(figura 1). Los experimentos se efectuaron 
en áreas relativamente planas, cubiertas 
con matorrales 


aislados y suficientemente accesibles, sobre 


parcialmente perturbados, 


todo seguros, lo cual permitió disponer de 
información adecuada y de buena calidad. 


Métodos 


En cada sitio se registró la precipitación total o 
incidente, la velocidad y dirección del viento, 
temperatura ambiental y humedad relativa cada 
veinte minutos, en una estación climatológica 
clásica (Vaisala modelo WXT510). Con base 
en el análisis de la composición florística 
de la vegetación, muestreos tipo Gentry y 
la determinación de los Índices de Valor de 
Importancia realizado por Becerril-Piña et 
al. (2007) fueron elegidos dos árboles de las 
especies mezquite (P. laevigata) y huizache (4. 
farnesiana), por ser especies dominantes de las 
microcuencas antes señaladas; el tamaño de 
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Figura 1. Localización de los sitios experimentales. 
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los árboles instrumentados se basó en la media 
dasométrica de los sitios de muestreo, en tanto 
que su ubicación se basó en que estuviesen 
suficientemente aislados, a fin de reducir los 
errores asociados con el efecto pantalla o sombra 
de árboles cercanos, y con la intención de 
observar la máxima capacidad de intercepción. 
A continuación se determinaron las medidas 
dasométricas de los seis árboles seleccionados. 
En el cuadro 1 se resumen las variables de cada 
uno de los árboles. 

El monitoreo de las variables indicadas en 
las ecuaciones (2) y (3) fue llevado a cabo con 
apego en la metodología descrita por Guevara- 
Escobar et al. (2007). En el caso del flujo de 
percolación (T,), la superficie de captación 
trató de englobar la mayor superficie posible 
de la proyección horizontal de la copa. Para 
reducir los errores asociados con el tiempo 
de retención y desalojo del flujo de la lluvia 
percolada, fue inducido un drenaje casi 
instantáneo, armando los paneles de captación 
con una fuerte pendiente de drenado, mayor 
que 15%. De igual forma se minimizaron 
los errores por temperatura, construyendo 
la superficie de captación con láminas de 
policarbonato transparente. La lluvia captada 
se almacenó en un depósito de plástico de 
250 litros, situado a un costado del tronco 
(figura 2). El flujo cortical S , fue monitoreado 
mediante la fijación de una canaleta alrededor 
del tronco o fuste, construida con una 
manguera de plástico, abierta por la mitad 
para eficientar la captación y el drenaje. El 
flujo captado se condujo hacia un garrafón de 
25 litros dispuesto en el piso. Es importante 
señalar que por razones económicas y técnicas 
fue difícil cubrir la forma irregular de la copa 
de los árboles para lograr una relación entre 
las áreas de 1:1. En el cuadro 1 pueden verse 
las proporciones de las áreas cubiertas con los 
paneles de acrílico. Gráficamente, en la figura 
2 puede observarse la distribución de las áreas 
formadas con los paneles respecto a la forma 
irregular de la proyección de la copa de los 
árboles. Ambos flujos, cortical y traslocación, 
fueron medidos al término de cada evento, 


transformado el volumen recuperado a lámina 
equivalente de lluvia. 

Las desviaciones de P,, originadas por la 
diferencia entre las superficies de captación 
formada de los paneles de acrílico y el área 
irregular de la proyección horizontal de 
la copa del 
suponiendo una distribución uniforme de la 
precipitación incidente. Con el supuesto de 


árbol, fueron corregidas 


un comportamiento lineal de los flujos, éstos 
fueron afectados con el coeficiente de superficie 
generado por el cociente del área del panel de 
acrílico y el área de la proyección en horizontal 
de la copa del árbol, como lo indica la expresión 
siguiente: 


pont =Spr Jen (10) 


donde A, es el área de la superficie de captación 
de los paneles de acrílico y A, es el área en 
proyección horizontal de la copa del árbol o 
sobre la superficie del suelo. Una vez corregida 
la precipitación neta fueron calculadas las E, 
aplicando la ecuación (3). En el cuadro 1 se 
enumeran por especie y por sitio los totales de P,, 
T, y S, correspondientes a los eventos utilizados 
en el cálculo de las E: 15 de El Carmen, 13 de 
Cadereyta y únicamente 10 de Amazcala, ya que 
por razones técnicas, en algunos eventos no se 
dispuso de las mediciones de T, y S,. 

Para la modelación de las E, fueron 
deducidos los parámetros requeridos por los 
modelos Rutter et al. (1975) y de Gash (1979). 
Para cada árbol se determinó la capacidad de 
almacenamiento (5) y coeficiente de percolación 
libre (p), aplicando el método de Leyton et al. 
(1967) y Rutter et al. (1971), respectivamente, 
este último con lluvias P¿<10 mm ante la 
carencia de las lluvias menores a un milímetro, 
como propone el método. Los parámetros $, y p, 
fueron deducidos empleando una metodología 
semejante a la establecida para el cálculo de 
la percolación libre (p) mediante la regresión 
5, versus P¿ Los parámetros D, y b del modelo 
de Rutter fueron derivados con la metodología 
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Figura 2. Esquema general de la instrumentación de los árboles y vista en planta de la proyección 
de la copa de cada árbol y del área de captación de T,. 
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de Schellekens et al. (1999). La evaporación 
potencial implicada en ambos modelos se 
calculó con la ecuación de Penman-Monteith, 
propuesta por la FAO (Allen et al., 1998). En el 
cuadro 1 se muestran los parámetros por cada 
modelo, en las columnas 12 a 20. 

Adicionalmente se exploró el efecto de otras 
variables climáticas en las E, mediante modelos 
de regresión lineal múltiple, explotando 
aquellas variables con el mayor grado de 
correlación, a fin de confrontarlo con los 
modelos de base física. Al mismo tiempo, fue 
analizada la robustez de la regresión por medio 
del criterio calibración-validación cruzada 
explorada por bloques. 

La estrategia de la modelación de las E, estuvo 
condicionada por el número de observaciones en 
los tres sitios. En estos términos, la calibración de 
los modelos de Rutter et al. (1975) y Gash (1979) 
se estableció de la siguiente forma: 


a) Calibración-validación cruzada por bloques. 

b) Medición de la bondad de los modelos para 
el periodo experimental con los estadísticos 
eficiencia, sesgo y la raíz del error cuadrático 
medio. 

c) Modelación de las E,para 2006 en Cadereyta 
y El Carmen, conservando los mismos 
parámetros utilizados para el periodo 
experimental. 


La bondad de los tres modelos, Rutter et al. 
(1975), Gash (1979) y la regresión lineal múltiple, 
se evaluó con los indicadores estadísticos 
siguientes: 


Eficiencia: 


1=n 


2 ( rl el q 


E=1-= (11) 


2 (Toy 


pos a 


Y a (12) 


Raíz cuadrada del error cuadrático medio: 


RMSE= q (y; model a (13) 


Resultados 


El 2006, año en el cual se realizó el estudio, 
puede considerarse como un año húmedo, ya 
que la precipitación en los sitios de Cadereyta 
(770 mm) y El Carmen (732 mm) se aproximó a 
los 800 mm, límite superior de la media anual 
reportada para la región de estudio. En El 
Carmen se registraron 54 eventos de lluvias y 
63 en Cadereyta. Durante los cinco meses de 
observación, agosto-diciembre, en Amazcala 
ocurrieron 24 eventos, produciendo 451 mm, 
aproximadamente el 60% de lluvia registrada 
en los otros dos sitios. En la figura 3 (a y b) 
se muestran las correlaciones por sitio y por 
especie (seis árboles) del tipo de la ecuación 
(8), E, versus P, correspondiente al comporta- 
miento a escala local. Individualmente, las 
regresiones resultaron con valores entre 
0.26<R?<0.87, siendo el huizache de Cadereyta 
el que presentó la correlación más baja. El rango 
promedio, las E, por evento, interpretadas a 
partir de la pendiente de la regresión, fue de 
10-38% en el huizache, y 12-20% en el mezquite. 
El flujo cortical resultó entre 1.9)<55.07%. En 
orden de magnitud, el de mayor contribución 
fue el huizache de Amazcala (6%) y el de 
menor, el mezquite de Cadereyta (1.9%). Es 
importante destacar que no obstante que el 
flujo cortical aparentemente es pequeño, es 
de vital importancia para el desarrollo de la 
vegetación, porque es agua disponible para la 
zona de raíces. 

El análisis ANOVA de las líneas de tendencia 
de las regresiones de P¿ versus E, mostradas 
en la figura 3 (a y b), demostró que no existen 
diferencias estadísticamente significativas 
entre ellas (P>0.05), lo que indica que los seis 
árboles presentan un comportamiento similar 
en cuanto a la cantidad de agua que interceptan, 
a pesar de la gran distancia entre sitios, superior 
a70 km. Debido a lo anterior, se decidió plantear 
dos opciones: la primera tratando a todos los 
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di Huizache a) Mezquite b) 
-=- A Cadereyta  E¡-0.984 + (0.096 * Pg) R2=0.26 === A Cadereyta E/1296+(0.199*P«) R2-0.87 
25] —— O ElCarmen Ej-0.119 + (0.196 * Pg) R2=0.78 — O ElCarmen E-1211 +(0.164* Po) R?059 
- m Amazcala  Ej=1.718+ (0.379 * P¿) R2=0.83 m Amazcala  E-0.541 + (0.121 * Pg) R20.44 
20 7 ñ 


Huizache c) Mezquite d) 
A Cadereyta A Cadereyta A 

25] O ElCarmen E;=0.740 + (0.251 * Pg) R?=0.62 O ElCarmen  EF1.398 + (0.146 * Pg) R*=0.50 
E  Amazcala EH Amazcala 


E, (mm) 


Pc (mm) Pc (mm) 
30 
Regional e) 
E¡=0.329 + (0.199 * Pg) R?=0.55 
25 y 
0] 
20 4 


Figura 3. Regresión entre la precipitación total y las pérdidas observadas por sitio, por especie 


y por unidad vegetal. 
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árboles de la misma especie como un solo 
individuo y la segunda tratando los seis árboles 
como unidad vegetal. 

La magnitud de las E, con base en el 
análisis por evento, fue del 15% para mezquite 
y del 25% para huizache. En la figura 3 (c y 
d) se presentan las correlaciones por especie, 
considerandolosseis árboles. La diferencia entre 
los promedios puede deberse a desigualdades 
morfológicas de ambas especies, ya que el 
huizache presenta un follaje más denso que el 
mezquite, permitiéndole interceptar más agua 
por unidad de superficie. 

En la figura 3e se muestra una correlación 
para todos los árboles integrados en una sola 
unidad vegetal, en correspondencia a una 
representación regional. La ecuación obtenida 
nos indica que las E, promedio de vegetaciones 
donde dominen estas dos especies son de 
20%. 

Para mejorar las predicciones de las 
regresiones lineares simples se analizó un 
modelo de regresión lineal múltiple (RLM) 
construido a partir de las tres variables con 
más influencia en las E: la precipitación total 
(P.), la velocidad del viento (u_) y la intensidad 
de la precipitación (); esta última modificada 
a la intensidad máxima del evento (imáx), en 
un lapso de veinte minutos, por ser el tiempo 
en que se presenta la mayor precipitación, 
comportamiento típico de las lluvias de 
la región. En el cuadro 2 se presentan los 
parámetros de 5 RLM a partir de diferentes 
combinaciones por bloques, y un resumen del 
análisis de variancia entre las E, calculadas por 
las 5 REM de cada árbol, de donde el huizache 
y el mezquite de Cadereyta son los únicos en 
presentar diferencias estadísticas entre los 
modelos. El análisis de sensibilidad de los 
modelos REM mostró que tienden a subestimar 
las E, en promedio, -0.5 mm, siendo menos 
eficientes para el huizache de Cadereyta y el 
mezquite de El Carmen. 

Determinada la sensibilidad de los cinco 
modelos RLM, se empleó el modelo construido 
a partir de todas las observaciones, para ser 
comparado con los modelos de Rutter y Gash, 


y modelar las E, para el 2006, ya que como 
se observó, el modelo RLM no presenta gran 
diferencia. 

Los modelos de base física, Rutter y Gash, 
fueron corridos con los parámetros indicados 
en el cuadro 1. Las E, simuladas en el periodo 
experimental y en 2006 con ambos modelos 
pueden observarse en la figura 4. En las 
columnas 3, 4 y 5 del cuadro 4 se indican los 
porcentajes de las E, en proporción con la 
precipitación total, correspondientes al periodo 
de observación y de 2006. El modelo RLM, 
sistemáticamente produjo mejores resultados, 
en comparación con los otros dos modelos, 
salvo para el huizache de Cadereyta. En las 
columnas de 5 a 13 del cuadro 3 se presentan 
los indicadores estadísticos de los tres modelos; 
específicamente, el modelo de Gash es el menos 
consistente, con valores de -4.2<E, <0.3. 

En la figura 5 se muestra la E, acumulada 
para el 2006, determinada con los tres modelos 
para las dos especies y los tres sitios; los 
modelos reproducen en forma semejante las E, 
acumuladas a lo largo del año. En las columnas 
6, 7 y 8 en el cuadro 4 puede observarse la 
comparación de los totales de los tres sitios y los 
porcentajes en proporción con la precipitación 
total. Con los modelos de Gash y Rutter, los 
rangos anuales varían entre 20-41% y de 17-29% 
con el modelo lineal múltiple. Sobresalen de los 
otros dos modelos las E, simuladas con Gash 
para el mezquite y huizache de El Carmen, 40 y 
41%, respectivamente. En el contexto regional, 
las E, estimadas con la regresión, mostrada 
en la figura 3e, resultaron de 23% (166 mm), 
22% (171 mm) y 22% (98 mm) para Cadereyta, 
El Carmen y Amazcala, respectivamente, las 
cuales corresponden al promedio de las E, 
representadas por especie (figura 3 c y d), y 
dentro del orden, las E, simuladas con los tres 
modelos. 


Discusión 
Precisando que los modelos de Rutter et al. 


(1975) y Gash (1979) tienen su origen en bosques 
templados con lluvias de larga duración, 
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Cuadro 3. Eficiencia (E), Sesgo (B) y raíz del error cuadrático medio (RMSE) para los modelos 
de Rutter, Gash y lineal múltiple. 


Rutter Gash Lineal múltiple 
Sitio Especie 
EW BY RMSES EW BY RMSE9 EW BY RMSE9 

(Dm (2) (5) (6) (7) (8) (9) (10) (11) (12) 13) 
Huizache -0.41 1.26 3.13 -0.09 2.32 2.99 -1.19 0.00 1.43 

Cadereyta 
Mezquite 0.57 0.44 2.46 -0.34 0.48 1.86 0.89 0.01 0.62 
Huizache 0.80 0.21 1.83 0.34 1.52 2.85 0.85 | 0.01 1.20 

El Carmen 
Mezquite 0.34 1.20 5.24 0.33 1.73 2.60 0.55 -0.89 2.79 
Huizache 0.30 8.25 9.19 En 1.34 4.01 0.78 0.02 1.83 

Amazcala 
Mezquite 0.39 3.24 3.58 -0.25 2.42 3.80 0.94 0.01 0.71 


WEficiencia, * sesgo, * raíz cuadrada del error cuadrático medio. 


Cuadro 4. Resumen de la modelación de las pérdidas por intercepción totales para el periodo 


Nota: () porcentaje respecto a la precipitación total. El número en negritas corresponde al porcentaje que se aproxima al promedio regional. 


experimental y el año 2006. 
Periodo experimental 2006 
Sitio Especie | Observados Gash Rutter RLM Gash Rutter RLM 
(mm) (mm) (mm) (mm) (mm) (mm) (mm) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) 
Hulradh 85 66 52 36 236 ie 131 
uizache 
(15) (29) (23) (16) (31) (16) (17) 
Cadereyta 
M ls 63 70 70 64 288 169 222 
ezquite 
q (27) (31) (31) (28) (37) (22) (29) 
En vil 94 74 72 297 161 169 
uizache 
(0) (27) (21) (20) (41) (22) (23) 
El Carmen 
M ls 75 102 94 76 291 190 177 
ezquite 
S id (22) (29) (27) (22) (40) (26) (24) 
> 
2 hiata 41 49 90 41 161 164 173 
Zi 
S (27) (33) (61) (27) (36) (36) (38) 
E Amazcala 
2 M e 25 40 45 25 148 88 113 
2 ezquite 
9 ii (17) (27) (30) (17) (33) (20) (25) 
Eo 
BS 
= 
3 


« vJB0/0u92], 


Mastachi-Loza, C.A. et al., Pérdidas por intercepción en mezquite (Prosopis laevigata) y huizache (Acacia farnesiana) de la región... 


Modelo (mm) 


Modelo (mm) 


Modelo (mm) 


Cadereyta Huizache 
O E¡Gash 
y E; Rutter 
MM E¡RLM 


4 


7 T T 7 Y pot r 7 r 7 7 Y 


2.0246 8 10 12 14 16 18 20 22 24 26 28 30 
Observación (mm) 


El Carmen Huizache 
O E Gash 
y Ej¡Rutter 
m E¡RIM 


2.0246 8 10 12 14 16 18 20 22 24 26 28 30 
Observación (mm) 


Amazcala Huizache 
¿0 E¡Gash 
y E; Rutter 
mM ERLM 


A A A A 
2.0246 8 10 12 14 16 18 20 22 24 26 28 30 
Observación (mm) 


Modelo (mm) 


Modelo (mm) 


Modelo (mm) 


Cadereyta Mezquite 


t 0 E¡Gash 
Y  E¡Rutter 
EM E¡RLM 


r 7 7 r 7 7 7 1 1 r y T 


4-2 0246 8 u10 12 14 16 18 20 22 24 26 28 30 


Observación (mm) 


+ El Carmen Mezquite 
+0 E¡Gash 

tw E;Rutter 

+ mm E¡RLM 

| a 
] ne E 


r 7 7 T T y T 7 7 T 1 T 7 T 7 


420246 8 1012 14 16 18 20 22 24 26 28 30 
Observación (mm) 


Amazcala Mezquite 
O E Gash 
wy E; Rutter 
mM E¡RLM 


420246 8 u10 12 14 16 18 20 22 24 26 28 30 
Observación (mm) 


Figura 4. Comparación de las pérdidas por intercepción, observación versus modelo Gash, 


Rutter y lineal múltiple. 
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Figura 5. Pérdidas acumuladas en 2006 en los tres sitios experimentales de huizache y mezquite, 
simuladas con los modelos de Rutter, Gash y RLM. 
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situación contrapuesta a las condiciones de la 
región de estudio, zona semiárida con régimen 
dominante de lluvias torrenciales de corta 
duración e intensas, se discuten los aspectos 
más relevantes de este trabajo. 

Adicionalmente, un sinnúmero de 
investigaciones señalan los errores atribuidos 
de la metodología utilizada en la medición 
de la precipitación neta, sin soslayar los 
efectos intrínsecos del medio observado: 
bosques cerrados o fragmentados, islas de 
vegetación y árboles aislados. La medición de 
las E, en árboles aislados suele realizarse con 
colectores aislados, líneas de canaletas y, en el 
mejor de los casos, con paneles. Sin embargo, 
como mencionan diversos investigadores, 
entre ellos, Crockford y Richardson (2000), 
es necesario una gran cantidad de colectores 
o canaletas para garantizar una medición 
correcta de E, para reducir el sesgo que 
genera la colocación aleatoria de colectores 
puntuales. Particularmente, este trabajo fue 
orientado a reducir los errores originados por 
la instrumentación en árboles aislados, con 
líneas de canaletas y colectores individuales 
(Chrockford y Richardson, 2000; Xiao et al., 
2000). 


empleada reduce la sobreestimación derivada 


Contextualmente, la metodología 
por las entradas laterales cuando la fuerza 
del viento incrementa el ángulo de la lluvia 
incidente (Guevara-Escobar et al., 2007). 

Entre los resultados a destacar del estudio 
está la semejanza con algunos trabajos de la 
zona semidesértica del norte México. Pese a 
que esta investigación fue realizada en la región 
semidesértica del centro de México, las E, (15- 
27%) en árboles aislados están dentro de las E, 
reportadas para individuos de la misma especie 
por Navar y Bryan (1994), Navar et al. (1999) y 
Carlyle-Moses (2004), cuyos estudios fueron 
enclavados en parches o islas de vegetación del 
semiárido del norte del país. Comparativamente, 
sus estudios reportan pérdidas por intercepción 
de 28% (Navar y Bryan, 1994), 18% (Navar et 
al., 1999) y entre 6-11% (Carlyle-Moses, 2004). 
Las discordancias entre los estudios pueden 
explicarse por la metodología de medición de 


do los patrones de la velocidad del viento, u; 
la intensidad de la lluvia, í; y la arquitectura de 
la vegetación. Sin embargo, refuerza la idea del 
modelo regional planteado, ya que éste propone 
unas E, del 20%, que caen dentro de la media 
observada de las E, reportadas en comunidades 
de matorral. 

Sorpresivamente, en esta investigación se 
encontraron valores similares de la capacidad 
de almacenamiento, S, con los valores repor- 
tados por Navar y Bryan (1994), entre 0.53- 
0.98 mm para huizache y 0.60-1.70 mm para 
mezquite. En la columna 12 del cuadro 1 se 
presentan las S de los seis árboles estudiados. 
La semejanza en S puede interpretarse como 
una expresión de una morfología similar 
entre las especies de A. farnesiana y P. laevigata 
del norte y el centro del país. Sin embargo, 
los valores del coeficiente de percolación 
en este estudio (p> 80%) en relación con los 
encontrados por dichos autores (23%) podrían 
indicar diferencias en el tipo de lluvia que 
se presenta en las dos regiones, producto de 
diferentes intensidades de precipitación. 

Regionalmente, la discrepancia entre los 
valores del parámetro de percolación, p, puede 
estar relacionada con la validez de método 
de Rutter et al. (1971). Dicho método se basa 
en observaciones de T, asociadas con lluvias 
ligeras, de poca magnitud y baja intensidad, 
para determinar el parámetro p. Al respecto, 
sólo los eventos utilizados en El Carmen 
presentaron las características mencionadas. 
Estas discrepancias pueden tener origen en 
una adaptación inadecuada de los modelos 
de Gash y Rutter para ambientes semiáridos; 
ante esta disyuntiva se planteó la opción de 
un modelo RLM. Este modelo consideró las 
variables de mayor nivel de importancia P,, u, 
e imáx; la reproducción de E, para los tres sitios 
de estudio, los seis ejemplares, los coeficientes 
de correlación resultaron relativamente altos. 
En comparación, la regresión lineal expresa 
que el conjunto de las variables u, e ¡máx 
afectan en un 16% las E, ya que localmente la 
velocidad ejerce una mayor influencia en las 
E, de árboles aislados, principalmente en su 
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distribución espacial en torno a la proyección 
horizontal del follaje, y la intensidad influye en 
el tiempo de saturación de la vegetación. Por 
otro lado, en promedio, la P, explica en un 55% 
la variabilidad de las E, logrando así, junto con 
u_ e imáx, explicar el 71% del proceso de las 
pérdidas por intercepción. 


Conclusiones 


El estudio se limitó a dos especies dominantes 
y los resultados pueden estimarse adecuados 
en el cálculo del balance hídrico de la región 
del semidesierto del centro del país y para 
reconsiderar su papel en la preservación 
de los ecosistemas semiáridos del centro de 
México. Los resultados obtenidos favorecen el 
conocimiento a escala global y regional de las 
pérdidas por intercepción; al mismo tiempo, 
podrían sustentar la calibración y validación 
de los componentes hidrológicos de los 
modelos de circulación global. A pesar de que 
las pérdidas por intercepción podrían jugar un 
rol menor en el balance hídrico de ecosistemas 
semiáridos y áridos, reduce la lámina de lluvia 
que genera el volumen de los escurrimientos 
superficiales, la infiltración o recarga profunda 
entre 15%< E,<30%, y el tronco de los árboles 
aporta del 3%<5,<5% del agua requerida para 
su desarrollo y preservación. Por otra parte, 
las E, regionales encontrados en este estudio, 
en promedio del 20%, ponen en evidencia 
la importancia que tiene la vegetación de 
la región semiárida del centro del país, 
específicamente el huizache (4. farnesiana) 
y el mezquite (P. Laevigata). A pesar de que 
ambos modelos, Rutter y Gash, sobreestimaron 
sistemáticamente E,, no presentaron diferencias 
estadísticas significativas entre las especies y 
los sitios, lo cual fortalece la hipótesis de base. 
De esta forma, el estudio llevado a escala local 
y por evento durante un solo periodo de lluvia 
constata una representación pertinente de la 
distribución regional (reducida variabilidad 
espacial) de las pérdidas por intercepción. Por 
otra parte, el modelo lineal múltiple evidencia 
que los modelos de base física deben considerar 


directamente los efectos del viento e intensidad 
de la precipitación en la determinación de la 
capacidad de almacenamiento y el flujo de 
percolación. 
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Abstract 


MASTACHI-LOZA, C.A., GONZÁLEZ-SOSA, E., BECERRIL-PIÑA, R. €; BRAUD, 1. 
Interception loss by mesquite (Prosopis laevigata) and huisache (Acacia farnesiana) in the 
semiarid region of central Mexico. Water Technology and Sciences, formerly Hydraulic 
Engineering in Mexico (in Spanish). Vol. I, no. 1, January-March, 2010, pp. 103-120. 


This study was performed in three sites of the semiarid region in central Mexico: “El Carmen” 
in the State of Guanajuato, and “Amazcala” and “Cadereyta” in the State of Queretaro. The 
main objective was to assess the interception loss process (E,) of the two dominant species 
in this region: mesquite (Prosopis laevigata) and huisache (Acacia farnesiana). For each 
site, two trees were selected, and the experiment was carried out following the methodology 
proposed by Guevara et al. (2007). Two physically-based models, Rutter et al. (1971) and 
Gash (1979), and one multiple regression model were used in order to model E. In 2006, 
precipitation in El Carmen and Cadereyta was 770 mm and 732 mm respectively, while in 
Amazcala it reached 451 mm from August to November. The interception loss observed for 
both species was very similar. The proportion of intercepted rainfall was 15%, 20% and 27 % 
for huisache and 27%, 22% and 17% for mesquite in Cadereyta, El Carmen and Amazcala 
respectively. The efficiency coefficient for the E, modelled during the study period for the three 
models was acceptable, ranging between -1.19<E,<0.94 and 0.62 mm <RMSE< 2.79 mm. At 
a region level, the E, was of 20%. 


Keywords: interception losses, evaporation, troughfall, stem flow, water balance, semiarid 
ecosystems, modelling. 
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